Resonating plaquette phase of a quantum six-vertex model 
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The simplest quantum generalization of the six-vertex model describes fluctuations of the order 
parameter of the d-density wave (DDW), believed to compete with superconductivity in the high- 
T c superconductors. The ground state of this model undergoes a first order transition from the 
DDW phase to a resonating plaquette phase as the quantum fluctuations are increased, which is 
explored with the help of quantum Monte Carlo simulations and analytic considerations involving 
the n-vector (n — 2) model with cubic anisotropy. In addition to finding a new quantum state, we 
show that the DDW is robust against a class of quantum fluctuations of its order parameter. The 
inferred finite temperature phase diagram contains unsuspected multicritical points. 

PACS numbers: 74.20.-z, 71.10.Hf, 74.10. +v 



The vertex models have unusual building blocks: ar- 
rows joined at a site forming the vertex. Nonetheless 
the statistical mechanics of these models are well de- 
fined, elegant, and have unexpected connections to more 
intuitively familiar models. For example, the transfer 
matrix of the classical six-vertex model is the XXZ- 
Hamiltonian [l| of Heisenberg spins and the correspond- 
ing transfer matrix of the eight vertex model is the XYZ- 
Hamiltonian. [2| Both of these models are solved by 
powerful mathematical methods with far reaching conse- 
quences. These vertex models arc not simply products of 
mathematical imagination but were originally proposed 
to describe phases of ferro and antiferroelectric materi- 
als. Q Thus they are effective descriptions of complex 
organization of matter. These vertex models are classical 
because they are made out of commuting variables. 

The interest in the quantum six-vertex model, to be de- 
fined below, is rooted in the unusual phenomenology of 
the cuprate high temperature superconductors. Q It has 
been argued that the origin of the pseudogap is an uncon- 
ventional broken symmetry in which a particle and a hole 
is bound in an angular momentum I = 2 state, result- 
ing in an order parameter, the d-density wave (DDW), 
which is effectively hidden. [j| The statistical mechanical 
description of DDW, which includes both thermal and 
quantum fluctuations of the directions of the bond cur- 
rents, is the quantum six-vertex model. |J] In particular, 
the description of DDW in terms of the six- vertex model 
resolves the puzzle as to why there are no associated spe- 
cific heat anomalies, as the phase transition corresponds 
to an essential singularity in the free energy. In this Let- 
ter we consider commensurate DDW as possible incom- 
mensuration wave vectors tend to be small in extended 
Hubbard models. 

Among the observed broken symmetries in the particle- 
particle channel are s, p, and ci-wave superconductors, 
where the Cooper pairs bind in the angular momentum 
channels / = 0,1, and 2 respectively. It is remarkable 
however, that in the particle-hole channel, the observed 



broken symmetries are mainly confined to s-wave sym- 
metry: s-wave singlet and triplet density waves, known 
as the charge and spin density waves respectively. Q 



Here we investigate, for the first time, the phase dia- 
gram of the quantum six- vertex model, which is an inter- 
esting and unusual statistical mechanical model in itself. 
A more practical goal is to examine the stability of DDW 
with respect to order parameter fluctuations. Until now, 
DDW has been treated within mean field theory; see, 
however, Ref. A full treatment of fermionic excita- 
tions is complex, but quenched fermions can be incor- 
porated as sources and sinks by enlarging the model to 
the eight vertex model. We show that DDW is remark- 
ably stable against quantum fluctuations. In addition, we 
note an unsuspected connection between the quantum 
six- vertex model and the (2 + l)-dimensional n-vector 
(n = 2) model with cubic anisotropy, which vindicates 
the numerical result that quantum fluctuations drive the 
DDW to a new phase, termed the resonating plaquette 
phase, via a first order transition. 



Our quantum six- vertex model bears some resemblance 
to the widely discussed quantum dimer model, but it does 
not exhibit the Rokhsar-Kivelson (RK) point 0, E| and 
differs from another recent attempt to quantize vertex 
models. 01 The latter, in which the RK point was built 
in by construction, had very different goals: the ground 
state degeneracy was exploited to provide, in part, an 
understan ding of topologically protected quantum com- 
putation. [lll|l2| Although our model contains topolog- 
ical sectors, distinguished by their total arrow polariza- 
tions in x and y- directions, the lowest energy belongs 
to the trivial topological sector of DDW, with zero total 
polarization (on an even-even lattice), see Fig. ^ For 
other interesting models involving constrained quantum 
dynamics, see Refs. 0,0]. 
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FIG. 1: (color online) Ground state energy per plaquette vs. FIG. 2: (color online) The squared DDW order parameter vs. 
t for different topological sectors labelled by the total arrow t for different lattice sizes. The inset is a blowup of the region 
polarization in the x and y-directions. (L = 12). near the phase transition. 



The Hamiltonian is 

h= Y, E{v)-t (®><0i + I0x0i) 

t?£ vertices p^plaq. 

(i) 

A configuration in this model consists of a set of inter- 
connected vertices that have four edges, arrows, which 
point either away from or into a vertex. In the six- 
vertex model only vertices with two arrows pointing 
into and two arrows pointing away from a vertex are 
allowed. The set of all possible vertex configurations 
forms the Hilbert space of the model. The first term 
in Eq. Q assigs a potential energy to a vertex con- 
figuration. We will choose E(*f*) = E(~^-) = and 

£(*f ) = E(~$+) = ) = £(4*) = !• The classical 
ground state is then DDW. The second term is quantum 
mechanical and reverses the arrows around elementary 
flippablc plaquettes . 

We investigate our model using a continuous-time dif- 
fusion Monte Carlo simulation yj| combined with the 
forward- walking technique ^| to extract ground state ex- 
pectation values. The accuracy of the simulations were 
improved by including a guiding wave function propor- 
tional to a power, a, of the number of flippable plaquettes 
in the evolved state. The value of a was determined by 
minimizing the variance of the local energy^] in short 
trial runs prior to the actual simulations. Simulation 
results for small systems with linear size L — 4 were 
checked against exact diagonalization results, and the 
largest system studied had L = 20. 

The DDW order parameter is defined as Oddw = 

IZ^EpGpiaql- 1 )^^^' where P*(Py) is the integer 
a;(y)-coordinate of the plaquette center p. C p sums the 
arrows around p. An arrow pointing in the counter- 



clockwise (clockwise) direction gives a contribution +1(- 
1). Fig. [3 shows (0ddw) as a function of i for different 
linear lattice sizes L. With increasing L the curves ap- 
pear to develop into a discontinuity at t — t c ~ 3.47, 
indicating a phase transition. 

To investigate the nature of this phase transition 
we construct histograms of Oddw consistent with the 
forward-walking method. The appropriate bin to incre- 
ment is selected from the value of Oddw occuring at a 
long time At = 3 in the past. A weight that corrects 
for the population bias needed to control the simulation 
is then added to this bin. The histogram is finally nor- 
malized by the total sum of weights. For small t the 
histograms are peaked at large positive and negative val- 
ues that move towards zero as t is increased. However 
this motion is not continuous. For t close to t r a 3- 
peak structure indicating a coexistence of phases jlj| is 
visible for the largest system sizes, Fig. When t is 
tuned through t c the peaks do not move, rather their rel- 
ative weights change, indicating a discontinuous change 
in Oddw m the thermodynamic limit. That the tran- 
sition is first order is also indicated by measuring the 
differentiated ground state energy per plaquette with re- 
spect to t. This results in numerical curves that develop 
into a discontinuity at the transition as the system size 
is increased, see the inset of Fig. [3J 

What is the ground state for t larger than t c l To an- 
swer this question we measure the density of flippable 
plaquettes m — (l/£ 2 ) J2 P E P - F p is 1 if the plaquette p 
is flippable, and is otherwise. We find that rif drops fast 
from 1 and saturates to a value nf ~ 0.56 for large values 
of t independent of system size. Thus the large-i ground 
state contains a significant number of flippable plaque- 
ttes. Another interesting quantity is the staggered den- 
sity of flippable plaquettes n sf = {1/L 2 ) £ (-1)^+^ F p . 
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FIG. 3: (color online) Histograms of Oddw at L = 20 for 
different values of t close to the transition. The inset shows 
the numerically differentiated (finite-difference) ground state 
energy with respect to t for different system sizes. 
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FIG. 4: (color online) The density of staggered flippable pla- 
quettes squared vs. t for different system sizes. The fig- 
ure indicates the resonating plaquette phase. Clockwise and 
counterclockwise arrangements of arrows are resonating on 
the "fat" plaquettes. 



A plot of (n 2 f ) vs. t for different system sizes is shown 
in Fig. 0] The staggered density is below the tran- 
sition and finite above. Thus there is a preference for 
flippable plaquettes on one of the two sublattices. At 
high values of t, \n s {\ w 0.39. For t > t c circulating 
arrows around plaquettes on one of the sublattices "res- 
onate" , changing their circulation direction constantly. 
We term this a resonating plaquette phase. This phase 
occurs also in a quantum six-vertex model with a built- 
in 

RK-point[3 

and in the Heisenberg model on a planar 
pyrochlore lattice. [2(| 

In order to identify proper order parameters for the 



quantum six-vertex model and infer the associated finite 
temperature phase diagram we write a six-vertex config- 
uration in terms of Ising spins at the center of the plaque- 
ttes. A vertex arrow points in the positive (negative) co- 
ordinate direction if its neighboring Ising spins are equal 
(opposite), see Fig.[S]upper right corner. By performing 
a 7r-rotation about the spin cc-axis for spins on the dashed 
diagonal lines in Fig. [S] the DDW shown can be identi- 
fied with the ferromagnetic state. In this new basis the 
Hamiltonian for the quantum six-vertex model becomes 
that of two inter-penetrating ferromagnetic Ising models, 
one on each sublattice, coupled by a frustrated four-spin 
term that couple Ising spins centered around a vertex. 
The couplings are chosen such that sources and sinks in 
the arrow-configurations are penalized energetically and 
a transverse field that flips the Ising spins plays the role 
of the quantum term. We replace the Ising spins with 
continuum fields Si, the label indicating the sublattice, 
and a term proportional to (Sf — l) 2 to preserve the spin 
magnitude. Thus, one arrives at the continuum action 
of an n-vector (n = 2) model with cubic anisotropy, ne- 
glecting irrelevant higher order derivatives : 



S 



0C 



dxo d x 




-Y 

2 ^ 



dSi 
dx 



(2) 



For our purposes, namely identifying the order param- 
eters, the precise values of the coefficients are not nec- 
essary. The isotropic n = 2 vector model with cubic 
anisotropy has been extensively studied. It has two 
ordered states for rn < 0. For v > it orders diago- 
nally. That is both Ising models (on both sublattices) 
have identical magnitudes of magnetization. For v < 
the n-vector model exhibits axis ordering: only one of 
the Ising models orders and the other remains disordered. 
The transition between these two ordered phases is first 
order. We identify the quantum phase transition in the 
quantum six-vertex model with this transition, and infer 
that the proper order parameters are the combinations 
of Ising order parameters on the two sublattices. 

We measure these order parameters in a simulation of 
Eq. n by assigning the value a = +1 to all Ising spins 
in the starting DDW-state. During the evolution a spin 
is flipped if the corresponding currents around the pla- 
quette is reversed, otherwise it is unchanged. The Ising 
order parameter on sublattice i is Ii = (2/ L 2 ) J2 P a v anc ^ 
takes values between —1 and 1. In Fig. [S] we have plotted 
I± = \I 2 ± if as functions of t. For t < t c i+ is finite 
and I— approaches zero for large system sizes indicat- 
ing identical magnetization magnitude on the two sub- 
lattices. For t > t c , both i+ and i_ approaches the same 
finite value (L — > oo), meaning that only one sublattice 
orders. The single-sublattice Ising order is an order by 
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FIG. 5: (color online) (/_} (pluses) and (7+) (circles) vs. t 
for different system sizes. Upper right corner: DDW configu- 
ration in terms of Ising spins. 



disorder phenomenon, as it is stabilized by the frequent 
plaquette flips, the resonating plaquettes, on the other 
sublattice. This makes the occurence of fhppable plaque- 
ttes on the ordered sublattice infrequent. 

As the temperature is increased, the action © will 
describe a two-dimensional n-vector model with cubic 
anisotropy, with renormalized values of u and v. For 
V = the transition is governed by an XY- fixed point, 
implying a Kosterlitz-Thouless (KT) transition above t c 
on the boundary between the ordered phases in the quan- 
tum six-vertex model. We will assume this also holds for 
small positive values of v and u in order to be consistent 
with the known fact that the t — thermal transition is 
KT-like.j^ It then follows by continuity that the ther- 
mal disordering of the DDW phase is in the KT univer- 
sality class for all < t < t c . For t > t c (v < 0) we also 
expect a KT transition in the vicinity of t c assuming that 
the renormalization group (RG) flows are still governed 
by the XY-point; however it is possible that the system 
will follow RG flows directed towards the fixed line ex- 
isting between the XY and the Cubic fixed point and 
that the KT-transition will change into a line of transi- 
tions with continuously varying exponents. On further 
increasing t this line might terminate and become first 
order consistent with the fluctuation-driven first order 
transitions known to exist in the n- vector model. [2l| We 
have sketched this possiblity in Fig. A direct transi- 
tion from KT-like behavior to first order transitions is 
also consistent with the RG flows for the n- vector model. 

The vertex models appear to be complex, but they are 
capable of encompassing unconventional order parame- 
ters such as orbital antiferromagnets, spin nematics, ex- 
citonic condensates, ice models, flux phases, and the pos- 
sible DDW state in the cuprates. Their quantum versions 
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FIG. 6: Possible phase diagram for the quantum six- vertex 
model as functions of flipping strength t and temperature T 
showing first order transitions (double lines), KT transitions 
(solid lines) and continuous transitions with varying expo- 
nents (dashed line). 



hold many surprises and are likely to exhibit fractional- 
ized excitations and non-trivial topological effects. In 
the present paper we have shown that they shed consid- 
erable light on the stability of the DDW against a class of 
quantum fluctuations and lead to a novel quantum state 
of matter, the resonating plaquette state. 

S. C. acknowledges a grant from the National Science 
Foundation: nsf-dmr 0411931 
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